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Abstract 

We present in details a numerical approach for solving supersymmetric 
quantum mechanical systems with a gauge symmetry valid in all fermionic 
sectors. The method uses a recursive algorithm to calculate matrix ele- 
ments of any gauge invariant operator in the Fock basis, in particular 
of the Hamiltonian operator, and can be used for any gauge group. We 
describe its application to a supersymmetric anharmonic oscillator model 
with discrete spectrum. 

1 Motivations 

For several years a revival of interest in Supersymmetric Yang-Mills Quantum 
Mechanics (SYMQM) can be observed. These supersymmetric quantum me- 
chanical systems can be obtained, most commonly, by a dimensional reduction 
of supersymmetric, D = d+1 dimensional, Af = 1 Yang-Mills quantum field the- 
ories to one point in space. Such procedure reduces the local gauge symmetry of 
the initial field theory to a global symmetry of the reduced quantum mechanical 
system. The physical Hilbert space of SYMQM is composed of states invariant 
under this global symmetry. As well all pertinent operators must be symmetry 
singlets. These constraints, being the remnants of the Gauss law, make the 
analytic construction of solutions nontrivial. 

The growing interest in these systems have several sources. On one hand, it 
is due to the their conjectured relation with a particular limit of M-theory [T], 
on the other hand, to the regularized dynamics of relativistic quantum mem- 
branes and supermembranes which they describe [2]. Moreover, their bosonic 
sector can be investigated as a zero- volume limit of 13 = 4, Yang-Mills quantum 
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field theory [3][1][S], providing results which can be compared with lattice cal- 
culations. Besides all this, the physically interesting features of SYMQM can 
be studied on their own. For these reasons much effort has been devoted to the 
evaluation of spectra of SYMQM. 

Among the variety of approaches, many numerical methods have been adapted 
to investigate supersymmetric Yang-Mills quantum mechanics and recently pro- 
vided new results. Some of them use path integral picture of quantum mechanics 
and Monte Carlo integration [6] [7] , other the Hamiltonian formulation of quan- 
tum mechanics and the Fock space methods [8] [9] [10]. Particularly, basics of 
the algorithm presented in this paper were already described in [IT][T2]. Our 
approach exploits the Fock space formulation, and therefore provides a non- 
perturbative way to calculate the eigenenergies and eigenstates. It was already 
applied to many systems uncovering interesting physics |13) . Although early at- 
tempts based on this principle proved to be very helpful in obtaining qualitative 
results, their potential was limited due to the rapid growth of the Fock basis. 
Comparing to them, the recursive algorithm which we present in this article, 
not only enables one to obtain numerical results up to very high precision, and 
thus, permits more quantitative considerations, but also makes the calculations 
in fermionic sectors possible. The main motivation for it comes from the study 
of D = 2, SYMQM. However, the method is much more general and provides a 
tool for evaluation of matrix elements of any invariant operator for any gauge 
group and in any dimension. Due to these multiple possibilities of extensions, 
this paper is the first of a series of articles presenting the results of studies 
of the SYMQM systems obtained with our numerical approach. Being the in- 
troductory paper, it contains, apart of the detailed discussion of the recursive 
algorithm itself, the summary of the whole framework which will be needed for 
the future work. Hence, in the following section we start with an introduction 
of the basic notions of SYMQM. Next, we construct the Fock space, which is 
necessary for the numerical calculations, and discuss some of its properties. The 
main part of the paper, the description of the algorithm, is divided into several 
sections. First of all, the calculation of the matrix of the scalar products is 
presented. Then, the orthonormalization procedure and the automatic evalua- 
tion of commutators and anticommutators is described. Only at that point the 
formula used for calculation of matrix elements of any operator can be clearly 
introduced. Eventually, the full recursive relations will be discussed. We finish 
by presenting a simple application of our method to a supersymmetric model 
with discrete spectrum. Conclusions and an outlook of future research directions 
will be provided in the last section. 

2 The framework 

In this section we describe several basic concepts constituting our framework. 
We start by introducing quantum mechanics in the cut Fock basis. Then, we de- 
fine the supersymmetric Yang-Mills quantum mechanics and derive a particular 
Hamiltonian which is studied numerically as an illustration of the approach at 



2 



the end of this article. In the remaining subsections we discuss the construction 
of the basis of SYMQM and its properties. 

2.1 Quantum mechanics in a cut Fock basis 

Quantum systems, which are described by a Hamiltonian operator expressed in 
terms of position and momentum operators, can be analyzed numerically in an 
eigenbasis of occupation number operators - a Fock basis [8'. Any occupation 
number operator can be written as a^a, where and a are bosonic creation 
and annihilation operators respectively, fulfilling the well-known commutation 
relations 

[flp, aq] = [al, aj] = 0, [op, aj] = Spq. (1) 

p and g in (H)) are indices which label the bosonic degrees of freedom. In the 
case of supersymmetric system we introduce fcrmionic creation and annihilation 
operators, / and respecting the anticommutation relations 

= {4,41 = 0, {fpjl} = 6pq, (2) 

where again p and q are indices which describe the fermionic degrees of freedom. 
Obviously, we also have 

[fp, aq] = [f;,aq] = 0, [fp, al] = [4, at] ^ q. (3) 

The momentum and position operators are expressed by creation and annihila- 
tion operators in the usual way 

Xq^ ^{aq + al), pq^ -^{aq~al), (4) 

enabling to express the Hamiltonian operator in terms of bosonic and fcrmionic 
creation and annihilation operators only. 

The construction of the basis starts with the definition of the Fock vacuum 
|0), as the state fulfilling the conditions 

aq\0) = 0, fq\0) = 0,yq. (5) 

Any other basis state can be obtained from |0) by a successful action of creation 
operators. 

Eventually, the action of the Hamiltonian operator, which is now an operator 
function of ^ and is straightforward in such basis. There is no conceptual 
difficulties in evaluation of its matrix elements, however such calculations may 
turn to be computationally demanding. The recursive algorithm described in 
this paper may be a solution to this problem. Once the Hamiltonian matrix is 
obtained, its eigenvalues correspond simply to the eigenenergies of the quantum 
system, and its eigenvectors to the eigenstates. 

The numerical analysis requires one last step, namely an introduction of 
a cut-off Ncut on the countably infinite Fock basis, since it is impossible to 
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deal with infinite matrices on a computer. There are many ways to introduce 
such a cut-off depending on the symmetries of the system being investigated. 
FinaUy, we have to perform calculations with several increasing N^ut and extract 
the physical results from the limit of infinite cut-off. The properties of such a 
procedure were analyzed in [S] [2] [H] [IS] . 



2.2 Supersymmetric Yang-Mills Quantum Mechanics 

The recursive algorithm will be described in the context of a relatively simple, 
D = d -\- 1 = 2 supersymmetric Yang-Mills quantum mechanics. In this way 
we will be able to focus our attention directly on the main features of the algo- 
rithm. However, it can be used for many different systems and it is conceptually 
straightforward to generalize it to the physically more interesting systems like 
D — A quantum mechanics [5]|12). In order to obtain the particular Hamiltonian 
which will be studied in the last section of this article, we will now derive the 
general Hamiltonian of I? = 2, SYMQM. To this end we extend the construc- 
tion of supersymmetric quantum mechanics |17j via the generalized creation and 
annihilation operators. We define the latter as 

A = + (6) 

At ^ l^{^p + w), (7) 

where the superpotential W = W{x) is a function of the position operator x. 
The supercharges are obtained from A and A"^ as Q ^ Af and — A^ . 
With the idea of SYMQM in mind, a simple generalization is to include a global 
SU{N) symmetry. Therefore, we postulate the supercharges to be given bj0 

Q= 7^ E {~'Pa+Wa)fa, (8) 
Q^ = ^ E (*Pa + W^a)/], (9) 



V2 



a=l 



where a is a color index of the adjoint representation of the SU (N) group, 
a = 1,...,N^ — 1, fa, fa are fermionic operators. Our system contains now 
N'^ — 1 bosonic degrees of freedom, described by Xa and pa , and A''^ — 1 fermionic 
degrees of freedom, described by fl and fa- The reduced Gauss law restricts 
the physical Hilbert space to only those states which are invariant under the 
SU{N) group. In order to facilitate the construction of such singlets it is useful 
to introduce a matrix notation. Thanks to the latter all singlets can be written 
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in terms of traces [H] . Therefore, we define 

a—1 a—1 

fi,j ^ 5Z -^a-^ij' -^'J ~ f'^'^iJ' 

a—1 a—1 

where Tf^ are the generators of the SU{N) group in the fundamental repre- 
sentation, i,j = 1, ... ,7V. Hence aU operators become operator valued matri- 
ces. We also introduce a simplified notation for a trace of any matrix, namely, 
tr(0) = (O). However, we will use this notation only when many traces occur 
and no confusion is induced. 

The Hamiltonian of such a system is given by the anticommutator of the 
supercharges, eqs. ^ and ©, and its general form read^ 

H = {Q\Q} = ^{paPa+WaWa) + ^idaWt + d,Wa)[fl fb]- (10) 

In the simplest case, we choose Wa = Xa and obtain a set of N'^ — 1 super- 
symmetric harmonic oscillators. For a slightly more complicated case, let us 
consider 

Wa = -dabcXbXc, (H) 

where g is the coupling constant and dabc is the totally symmetric tensor of the 
SU{N) group. The Hamiltonian ea.([TU)) reduces to, 

H ]^{PaPa + ^dabedecdXaXbXcXd) + l^dabcXaUl Jc]- (12) 

One of the interesting features of this model is a nontrivial bosonic potential of 
fourth order which for the SU (3) group is simply 



dabedecdXaXbXcXd^\^{^^{Xaf^ ■ (13) 



Therefore the system of eq. (|T2|) is expected to have a discrete spectrum of bound 
states. In order to rewrite this Hamiltonian in terms of traces we use the 
following identities valid for any SU (N) group, 

1 2 
^PaPa = tr p , 

~X(iXi)XQXddabedecd — tl X ~ (^^ ) ' 

\xafUcdabc = tr(a;/V) - tr(a;//t), 



We adopt the notation when a repeated index is assumed to be summed over. 
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and eventually obtain 

H^trp^+ g^{tY - ^(tr x^f) + 2g tr(x[/t, /]). (14) 

We suspend the detailed discussion of this Hamiltonian to section |4] where the 
numerical results obtained with the recursive algorithm will be presented to- 
gether with some analytic calculations. In the following subsection we will 
describe the construction of the Fock basis for SYMQM and some of its prop- 
erties. 

2.3 Elementary bricks and Fock basis of the SU{N) 
SYMQM 

The Fock states are eigenstates of some occupation number operators. In the 
case of SYMQM models we consider gauge-invariant occupation number oper- 
ators, 

tr a^a = ^ at'a', (15) 

q=l 

Most of the Hamiltonians that we have considered so far conserve the fermionic 
occupation number. Therefore, it is physically motivated to construct the Fock 
basis independently in each subspace of the physical Hilbert space with a def- 
inite fermionic occupation number. Moreover, states containing different total 
number of quanta are orthogonal, so we can further divide the fermionic sec- 
tors into subspaces with given number of bosonic quanta. However, usually the 
Hamiltonian does not conserve the bosonic occupation number and thus mixes 
different bosonic subspaces. We start the construction of the Fock basis in the 
purely bosonic situation, and then turn to the fermionic sectors. 

2.3.1 Bosonic elementary bricks 

A general eigenstate of the occupation number operator, having ub quanta, can 
be written as [H] 

Ins) = E '?^^i.»......^..«---«ljO), (17) 

where T is a group invariant tensor. It can be shown |18j that any such invariant 
tensor can be expressed as linear combination of products of trace tensors. The 
latter are just traces of products of matrices, of which the simplest ones are 
tr T'^T'' = i(5afc and tr T^T^T" = irf^bc + jfabc- Therefore, in matrix notation 
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a state |ns) can be rewritten as 

Wb) = J2 7fc.,...,fc„,(at2)'-nat3)^-3...(at"-)'=".|0), (18) 

where Jk2....,k„^ are arbitrary coefficients and the sum is over aU such combina- 
tions of exponents kj that, X]j=2i% ~ ""S, so that the state |ns) is composed 
of ub quanta. Notice that even though a single quantum created by a| ^ is gauge 
dependent, total numbers of quanta, n^, being the eigenvalues of the operators 
eqs. p5t and (flH)) . are gauge independent. 

Subsequently, Ins) can be significantly simplified with the use of the Cayley- 
Hamilton theorem. It states, that any matrix. A, satisfies its own characteristic 
equation. Therefore, we have for the SU{2), SU{3) and SU{i) groups 

SU{2): 1(^2)1 = 0, 

SU{3): A^~^-{A^)A~^-{A^)I = 0, (19) 

SU{A) : A^ - 1(^2)^2 _ 1(^3)^ ^ 1 (^2)2^_ l(^4)j ^ 0^ 

One can use these equalities to reduce traces containing more than N operators 
of the same kind to simpler ones. We demonstrate this on an example with 
A — a'^ and some arbitrary operator B, which can be any operator involving 
bosonic or fermionic creation and annihilation operators. Particularly, B can 
be again a single bosonic creation operator. Thus, multiplying eqs. (fT9|) by B 
from the right-hand side and taking the trace, we obtain a set of relations, such 
as 

SUi2) (atflts) = ^{a^a^){B), 

SUiS) (atfltflt^) = -(atat)(at5) + -{a'< a'' a^){B), (20) 
2 3 

B) = ha'^a)){a)a'^B) + \{a) a^){a) B) - Ua'^a'^fiB) 
2 3 8 

+ ^{a^a^a^a^){B). (21) 
Hence, a general state with ng quanta for some given N ^ simplifies to 

\nB)N = J2 7fc.,...,fc„ (a^YM«^')'^---(«^^)'"|0), (22) 

{ Ef=2 J'=3="b} 

where the traces with more than N creation operators were reduced and the 
highest trace is now (a^^). 

We are now in position to define the set of bosonic elementary bricks, which 
is the set of A^ — 1 linearly independent single traces of creation operators, which 
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SU{2) 


SU{3) 


SU{4) 




(fltatat) 


(fltatfltat) 



Table 1: Elementary bosonic bricks for SU{2),SU{3) and S'C/(4). 



cannot be further reduced by the Cayley-Hamilton theorem. Table [T] contains 
examples of such sets for = 2, = 3 and = 4. Products of powers of 
elementary bosonic bricks acting on the Fock vacuum compose the set of states 

|(atYH«^Y^---(«^^)'"|0)}^„ ^ ^\{nB}) (23) 

which spans the subspace of the Hilbert space with ub bosonic quanta. We 
adopted a generalized notation in which llns}) is a vector of all states with ns 
quanta. The set of states eq. (P5|) is excessively called the Fock basis, though 
in general it is not orthonormal. Only after the application of an orthonor- 
malization procedure it will be transformed into a basis. By acting on it with 
additional elementary bricks and orthonormalizing, one can obtain the basis in 
sectors with yet higher number of bosonic quanta. In this way, starting with 
the Fock vacuum |0), one can recursively generate the basis for any ub- 



2.3.2 Fermionic bricks 

The definition of fermionic bricks is a bit more complicated. Besides the fermionic 
bricks which are single-trace operators, as, for example, all fermionic operators 
in the case of the SU{2) symmetry group (see table [2]), we must also take 
into account bricks which are multiple-trace operators. They appear in higher 
fermionic sectors in the case of symmetry groups with N > 2 (see tables |3] and 
1^ . In order to make this distinction clear we will now introduce some definitions 
and appropriate notation. 

In analogy to the set of elementary bosonic bricks, we define the set of el- 
ementary fermionic bricks. The latter will contain all single traces with np 
fermionic creation operators, which cannot be further reduced by the Cayley- 
Hamilton theorenjfl. An algorithmic way to obtain it, is to start with the set 
of bosonic elementary bricks for a given A^. Then, for each trace, one should 
perform np times either one of the following operations: replace one bosonic 
creation operator by a fermionic one or insert one fermionic creation opera- 
tors into the trace. In order to obtain a complete basis, the set of elementary 
bricks must be enlarged by operators, which are products of fermionic elemen- 
tary bricks with smaller number of fermionic quanta and contain np fermionic 

•^There is no simple counterpart of the Cayley-Hamilton theorem for anticommuting ma- 
trices. However, an appropriate choice of the B operator in H21|l will produce identities which 
can be used to simplify or exclude, due to linear dependence, some of the possible composite 
fermionic bricks |19| . 
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np = 1 


np = 2 


np = 3 


(r«t) 




iPPP) 



Table 2: SU{2) fermionic bricks. 



creation operators in total. In this way we ensure that all possible invariant con- 
tractions of np fermionic creation operators with a number of bosonic creation 
operators are considered. However, the problem of linear dependence appears 
and one has to pick out only the linearly independent operators. The linear 
independence can be checked by an explicit calculation of the determinant of 
the Gram matrix of states constructed with those bricks. Fortunately, there 
exists also an independent, and more direct way of computing the number of 
Fock basis states which will be described in subsection 12.41 

The enlarged set of all linearly independent fermionic bricks will be refered 
to as the set of composite fermionic bricks. Its elements will be labeled by an 
index a, and denoted by C^{nB,np,a), where ns and np are the number of 
bosonic and fermionic creation operators, respectively. Such notation is used in 
order to treat all bosonic and fermionic bricks in the same way. The additional 
index a distinguishes the operators with the same np and n^- 

Moreover, we will denote: 

• the number of composite bricks with n p fermionic and ns bosonic quanta 
by d{np,nB), 

• the total number of composite bricks with np fermions by d{np). 
Obviously, 

d{np) = y^d(n_F,nB). 

We extend our notation to the set of bosonic elementary bricks 

(at"B) = ct(ns,0). 

As an example. Tables 0, Q and (0]) contain the sets of composite bricks for 
N = 2, N = 3 and iV = 4 for some fermionic sectors. Let us briefly comment 
on table [H In this simple case, there is exactly one fermionic operator in each 
fermionic sector (see section r2.4p and they are just elementary fermionic bricks. 
One can prove, by simple arguments, that other operators are not possible. For 
example, an operator of the form (/^a^a^) vanishes identically, since it must 
involve the totaly antisymmetric tensor e'-'*' summed with a symmetric combi- 
nation of bosonic operators a'^a'^'^. In analogy, the possible operator (/^a^)^ 
also vanishes, since it is a square of an anticommuting operator. Such reason- 
ings facilitate the explicit construction of the set of composite fermionic bricks 
for gauge groups with N > 2. 

Once the set of composite fermionic bricks is constructed, it is easy to write 
down a general state with ns bosonic and np fermionic quanta for a given 
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np = 1 


np — 2 


np = 3 


rii? = 4 




(/tatflt/tat) 
(/tat)(/tatat) 


iPPPa^) 
(/^at)(/t/tat) 

(/tatat)(/t/tat) 
(/tatat)(/t/tatat) 


(/^a^at)(/t/t/t) 

iPa^)ia^PPP) 
(/V^at){/t/tat) 
(/tatat)(/t/t/tat) 

(/tat)(/tatat)(/t/tat) 
{PPa^){pa^pa^ar} 



Table 3: SU (3) fcrmionic bricks. 



gauge group SU{N). One has to take a linear combination of states obtained 
by applying one of the composite fermionic bricks with n p fcrmionic quanta to 
a general bosonic Fock state eg. ([2^ . Hence, 

d{nF ) 

\nB,np)N = ^ C\n,np,a)x 

a=l 

E 7fe.,...fc„(a) Ct(2,0)'=^Ct(3,0)'=^...Ct(iV,0)'="|0), (24) 

where the coefficients 7fe2....,fcN (o^) '^^^ depend now on a. In order to ensure 
that the total number of bosonic quanta is ns, we have apply the operator 
C'^ {n^np^a) containing n bosonic creation operators, to a purely bosonic state 
with ub — n quanta. In analogy to the bosonic case, we can define the set of 
states, 

{c\n,np,a)C^{2,0)'^ ...CHN,0)'^-\0)} =\{nB,np}) (25) 

which after orthonormalization will give the basis in the subspace of Hilbert 
space with ub and np bosonic and fermionic quanta, respectively. The linear 
independence and completeness of the set of composite fermionic bricks ensures 
that eq. ipS]) form indeed a complete set of states in the fermionic sectors. 

2.4 Counting the number of SU{N) SYMQM Fock basis 
states 

In the preceding section we have shown, that the Fock basis in any sector is 
given by eg. (1251) or eg. ()25l) . As an example, table [5] contains few simplest states 
of the Fock basis with the SU{2) symmetry. 
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np — 1 


np = 2 


(/tatatflt) 


(/^/^«^) 
(/t/tatat) 

(/t/tfltatat) 

(/tat/tfltat) 

(/tat)(/tatat) 

(/tat/tatatflt) 

(/tat)(/tatatat) 

(/tatat/tatatat) 

(/tatat)(/tatatat) 



Table 4: SU{4) fermionic bricks. 





np = 


nj? = 1 


nj? = 2 


rii? = 3 





|0) 






(/VVnio) 


1 




(/^«^)|0> 


(/T/tat)|0) 




2 


|2)^(«tat)|0) 






(/W^)|2) 


3 






(/W)|2) 














2n 


|2n) = (atat)"|0) 






(/VVn|2n) 


2n+l 




(/tat)|2n) 


(/t/tat)|2n) 















Table 5: Construction of the basis for the SU{2) gauge group. 



Before considering the set of vectors as in table (O as a basis one must check 
its hnear independence. This can be done by explicit calculations of the deter- 
minant of the Gram matrix however such computations become cumbersome for 
larger number of states. Fortunately there exist an alternative, group-theoretical 
way of calculating the total number of linearly independent gauge-invariant Fock 
states. It was suggested by Janik and elaborated by Trzetrzelewski |19) |20j. 
Instead of an explicit construction of basis states, this approach exploits the or- 
thogonality of the characters. Such an alternative method is of great practical 
value since it may serve as a crosscheck to our recursive algorithm. 

2.4.1 Character method 

Let D{nB,np) be the number of gauge-invariant, linearly independent states 
with ub bosonic and np fermionic quanta. D{nB,np) can be obtained from 
the orthogonality relation of the characters of the SU{N) group. 

Each bosonic and fermionic creation operator transforms according to the 
adjoint representation of the SU{N) group. Hence, the products of creation 
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operators, which are needed for the construction of basis states, transform as 
products of the adjoint representations. From the representation theory it is 
known, that the square of any irreducible representation is reducible and can 
be expressed as a sum of a symmetric and antisymmetric parts. This statement 
written in terms of characters reads, 

x(i?)xx(i?) - MR) >^ x{R)] + {x{R) >^ xiR)} 

= l{x'iR) + xiR')) + l{x'iR)-xiR')), (26) 

where the symbols [. ..] and {. • ■ } denote the symmetric and antisymmetric 
part, respectively, and x(i?^) is the trace of the matrix of the representation R 
squared. A generalization of the eq. ipSj) is known as the Frobenius theorem pT] . 
and has a practical meaning, since the characters of powers of R are often ex- 
plicitly given. It gives the expressions for the symmetrized and antisymmetrized 
characters of a product of p representations R, 

ixL.x(«)i = E rii^^. p^) 

where the sum is over all partitions of the number p into numbers 2, . . . ,p, ij 
being the multiplicity of the j number in a given partition. Thus, the most 
general product of ns bosonic and np fermionic creation operators will be in 
the representation, which character is equal to [><k=iX{R)]{^k=iX{R^)}^ where 
R denotes now the adjoint representation of the SU (N) group. From the or- 
thogonality property of the characters we have 

D{nB,nF) = J 1 [xlZiX{R)]{y<ZiXm, (29) 

where 1 stands for the character of the trivial representation and 'i/i5[/(Ar) is the 
group invariant measure on SU{N). 

A convenient parametrization of the group manifold is by N'^ — 1 Euler angles 
ai, all defined on [0, 2tt]. For example the group elements of SU{3) read 

JJ — giA3aigiA2a2giA3Q3giA5Q4giA3a5giA2Q6gi>'3a7gjA8Q8 i'iO) 

and the generalization of the above equation to SU{N) can be found in [53]. 
The last element needed to calculate D{nB,np) are the characters x{R^)- They 
are given by the Weyl formula [24] 

JV N 

X{R) = ^'^"'""^^ - 1' X(i?') = E e^'^"-""^^ - 1- (31) 
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The invariant measure reads |19] . 

^.-(«) - M - e-or^p(E"0> (32) 

i—l i<j i—1 

where Sp{x) is a periodic deha-function given by 

N oo N 

^P(E"0= E '5(E««-27rfc), (33) 

i—l k— — c>o i—l 

with k integer. 

2.4.2 Generating functions for D{nB,np) 



Eq. ([29|) is difficult to evaluate for any N, however it was calculated for few 
simplest groups [19 . In these cases, the numbers D(nB,np) can be encoded in 
a generating function with two parameters t and s, G{t, s), 

G{t,s)= E D{nB,nF)r-{-sr-, (34) 

nB,np 

which is very useful in practical applications. For iV = 3, G{t, s) can be ex- 
pressed in terms of simple polynomials in t 19], namely 

^(*' ^) = ( n T3^) E (35) 

fc=2 i=0 

and the polynomials Ci(t) read, 



co{t) 


= 1, 


Cl(i) 


= t + 


C2(i) 


= t + t^ + 2f^, 


C3(i) 


= 1 + i + 2<2 + 3t^ + i'' 


C4(i) 


= 2t + + + 2t'^ , 


Cs-i{t) 


= c,(t). 



(36) 

In this form some information contained in G(i, s) become evident. The term 
proportional to is equal to the generating function for the number of partitions 
into numbers {iV, iV — 1,...,2}. Obviously, there are as many states with 
ub quanta as there are ways of obtaining ub from multiples of the numbers 
of quanta contained in the elementary bosonic bricks. Furthermore, from the 
polynomials Ci (t) the combinatorial interpretation of the multiplicities of states 
in the fermionic sectors can be read off. Particularly, the number d{nF,nB) of 
composite fermionic bricks with ub bosonic quanta in a given fermionic sector 
is simply given by. 



1 d"^ 

d{np,nB) = ^_^c^At) 



(37) 

t=o 
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As an example, let us take the polynomial ci{t) for the SU{3) group. We have 
one brick with a single bosonic quantum - the (a^/^) brick, and one brick with 
two bosonic quanta - the (a^a^/^) brick. A less trivial example is given by the 
C2 (t) polynomial. Apart of the two operators with a single and double bosonic 
quanta, (a^f^f^) and {a^a^ f^), respectively, we now have two operators with 
three bosonic quanta, namely, {f^a^a^ f^a^) and {f^a^){f^a^a'^). 

Thus, we can crosscheck the number of basis states obtained by direct con- 
struction and elimination of linearly dependent states with the one computed 
using the above group-theoretical predictions. Table [6] presents the multiplicity 
of basis states with given and np quanta for iV = 3 up to ng ~ 30, cal- 
culated with both methods. We simply give a single set of numbers since the 
results agree exactly. 

2.5 Symmetries of the SU{N) SYMQM Fock basis 

The Fock basis of SYMQM, constructed in the way described above, has several 
nontrivial symmetries which can be interpreted as announcements of dynamical 
symmetries of the Hamiltonians of SYMQM. The most important of them are 
the supersymmetry and particle-hole symmetry. The former can be observed as 
a matching of eigenenergies from neighboring fermionic sectors, while the latter 
is defined as a matching of spectra from the sector with p fermions and the 
sector with iV^ — 1—p fermions with < p < N'^ — 1. Let us now describe three 
observations jT5] of nontrivial relations among the multiplicities D{nB,nF) and 
their interpretations in terms of these symmetries. 

• Supersymmetry 

— For each ub there are as many bosonic basis states (states with np 
even) as fermionic basis states (states with np odd). 

Vni, ^ D{nB,nF)^ ^ D{nB.nF)- 

np - even np - odd 

Its validity can be checked explicitly for SU{i) case, either in table 
[5] by summing the numbers of states with n p even and odd in each 
row separately, or in eqs. (fM|) by summing appropriate polynomials. 
This relation can be also exactly proved for any N using the general 
form of the generating function [19j . 

— Summing the number of states along diagonal lines with Ub + np or 
Ub — np fixed, for np even and np odd separately yields the same 
results, 

VriB ^ D{nB±np,np)^ ^ D{nB ± np,np), 

rip - even np - odd 

Such diagonal lines correspond to the action of the supersymmetric 
gauged harmonic oscillator supercharges, Q — (fa), and = (/^at). 
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Table 6: Multiplicity of Fock basis states with given number of bosonic and 
ferniionic quanta for the SU (3) gauge group. These numbers can be obtained 
from the generating function eq. p6p and from the orthonormalization proce- 
dure of the algorithm independently. 



They are of special interest since introducing the cut-offs in the con- 
secutive fermionic sectors according to one of these lines allows to 
obtain an exact supersymmetric degeneracy for finite cut-off. Again, 
this result can be proved for any N using the generating functions 

m- 

This is still not true supersymmetry. There is no dynamical supermulti- 
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plets, because at this stage we have not defined any Hamiltonian. However, 
it is interesting that ah'eady at this level we have such a matchings. 

• Particle-hole symmetry 

— For any n^, the number of states in the sector with np fermions is 
equal to the number of states in the sector with iV^ — \ — np fermions, 

V„^ D{nB.nF) = D{nB,N'' -l^np), np = Q, . . . , - I. 

If we expect that the spectrum in the sector with np fermionic quanta 
coincides with the spectrum in the sector with N"^ — 1 — np then the 
equality of the multiplicity of basis states in those sectors can be 
interpreted as a nontrivial announcing of the particle-hole symmetry. 

Summarizing, already at kinematic level one can find symptoms of the sym- 
metries of the SYMQM systems. 

3 Description of the algorithm 

The main idea of the recursive algorithm has been already described in [TT] [I^] . 
It relies on the observation that the most efficient way to evaluate a matrix 
element of an operator is to relate it to simpler matrix elements of some op- 
erators, which have been already evaluated at an earlier stage of calculations. 
In this way the explicit construction of the Fock basis vectors is not necessary. 
As an input the algorithm needs the commutators/anticommutators of elemen- 
tary bricks and any other invariant operators, which appear in these resulting 
commutators /anticommutators. 

In order to expose the algorithm in a clear way, we will start by explaining 
the construction of the Fock basis and then the calculation of the matrix of scalar 
products. Having such a matrix, one can orthonormalize the basis vectors. The 
procedure used to this end will be described in the subsequent subsections. At 
that point the formula used for calculation of matrix elements of any operator 
will become evident. Eventually, the full recursive relations will be presented. 

3.1 Recursive construction of the Fock basis 

Although the Fock states are not explicitly needed for the computation of matrix 
elements of an operator, they do appear in the labeling of those matrix elements. 
Hence, we should have a recursive way of obtaining the Fock basis. Let us 
assume that such basis is already constructed in the sectors with the number 
of bosonic and fermionic quanta smaller than ub and np respectively. Then, 
the Fock basis in the sector with ub bosonic and np fermionic quanta can be 
created as the sum of all states obtained by the action of appropriate bricks on 
the already generated Fock basis states. In our generalized notation this can be 
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written as 

N 

\{nB,nF}) = ^ C''{k,0)\{nB - k,nF})- 

k=2 

Note that in general suclr states will not form an orthonormal set of states. 
Moreover, the same state may appear in several copies, differing in the order of 
successive bricks used to build it. Those duplicates will be treated as distinct 
states. The basis is obtained once this redundancy is removed and the remaining 
states orthonormalized. 

3.2 Matrix of scalar products 

The Fock basis obtained recursively form a complete set of states which are not 
orthonormalized. Thus, one has to calculate the matrix of scalar products. It is 
sufficient to calculate the scalar products among the states containing a given 
number of quanta, and n p , since those having different number of bosonic 
or fermionic quanta are orthogonal by definition. The matrix of such scalar 
products will be denoted by S{nB, np)- 

The definition of S{nB,nF) diff'ers from the standard definition of the Gram 
matrix in few aspects. The traditional Gram matrix contains the scalar products 
of linearly independent states and for our systems it is a D{nB, np) x D{nB,np) 
matrix. Contrary, since in our algorithm the Fock basis is defined recursively, 
the S{nB,np) is the matrix of scalar products of states that have ub bosonic 
and np fermionic quanta and are obtained by the action of appropriate bricks 
on states from sectors with smaller number of bosonic and fermionic quanta. 
Therefore, the matrix S{nB,np) is usually bigger than the Gram matrix since 
some of the states can be included several times. Such redundancy is best 
illustrated by the following example. For the SU (3) gauge group, there is only 
one Fock state with 5 bosonic quanta, |{5, O}) = C'l'(2, 0)Ct(3, 0)|0). However, 
the calculation of the matrix of scalar products 5(5, 0) yields a 2 x 2 matrix 
instead of a single number, namely, 

s(5,o) = ({5,o}|(5,o})^(f(|2)»;> 

where the diff'erent matrix elements 5'(5,0)ij correspond to the multiple possi- 
bilities of pulling out an elementary brick out of the state |{5, 0}), 

|{5,0})=Ct(2,0)|a)=Ct(3,0)|/3), 

where |a) and are appropriate remaining states. Such doubling must be 
eliminated. The procedure which we used to achieve this will be described in 
the next subsection. 

The evaluation of S{nB,np) can be divided into two separate cases. 
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If np — 0, we pull out one bosonic elementary brick from each of the basis 
states. Since we have A^ — 1 bosonic elementary bricks (see tables [2 El E]), there 
is iV 1 different ways to do this, if only ub > N. If ub < N we can pull out 
only Ub — 1 different bosonic elementary bricks. Thus, in a generic situation, 
S{nB: 0) will be a (TV — 1) X (A^ — 1) matrix of the form. 



SinB,0) 



I S{nB, 0)2,2 S{nB, 0)2,3 ■■• 5(713,0)2,^-1 ^ 
S{nB, 0)3,2 S{nB, 0)3,3 ■•■ S{nB,0)3,N-i 



\ S{nB,0)N-l,2 S{nB,0)N-l,3 ■•■ S{nB,0)N-l,N-l J 



where 

S{nB,0)p,, = {{riB ~ p, 0} \C{p, 0)C^(g, 0)\{nB - g, O}). 

The matrix element, calculated by extracting the C{p, 0) elementary brick from 
the left state and the C''{q, 0) elementary brick from the right state, denoted by 
S{nB^ 0)p,q, can be expressed in terms of matrix elements of operators between 
basis states with lower number of bosonic quanta as, 

5(ns,0)p,, = {{nB-p,0}\[C{p,0),C\q,0)]\{nB-q,0}) 

+ {{nB~p,0}\C^{q,0)Cip,0)\{nB~q,0}). (38) 

Hence, we expressed S from the sector with n b quanta in terms of matrix ele- 
ments of operators evaluated in the sector with ub — q quanta. The procedure 
to calculate these matrix elements is described in subsection [ 



If Ui? 7^ 0, we first want to anticommute the composite fermionic bricks with 
Up fermionic quanta. Again, in a generic case, we treat S{nB,nF) as a matrix 
of size dlnp) x d{np) and of the form 

SinB.np) = 

I S{nB,np)i^i S{nB,np)i^2 •■• ^(tib, 71^)1, \ 

S{nB,np)2,i S{nB,np)2,2 ■■■ S{nB,np)2,d{nF) 

\ S{nB,np)j_(^nF),l S{nB,np)j_(^nF),2 •■• 'S'(nB, nj;-)rf(„p)^d(nF) / 

where 



S{nB,np)p,q = {{ub - np,0}\C{np,np,p)C^{nq,np,q)\{nB ~ nq,0}). 

Let us consider one of the scalar products, denoted by S{nB,np)p,q, and ob- 
tained from the scalar product of two basis vectors, first containing C{np, np^p), 
second containing C'' {nq,np, q). We have 

S{nB,np)p^q = {{ub - np,0}\{C{np,nF,p),C''{nq,np,q)}\{nB - nq,0}). 

(39) 
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The anticommutator {C{np,np,p),C^ {nq,nF,q)} is a normally ordered oper- 
ator containing only bosonic creation and annihilation operators. If it is not, 
one can always bring it to such forrrQ. The matrix elements of the operator 
C^{nq, Hp, Q)C{np, np,p) vanish since there are fermionic annihilation operators 
acting on the Fock vacuum. Thus, by anticommuting all fermionic operators 
in one step, we can express S{nB, np) in terms of matrix elements of operators 
between states from the bosonic sector exclusively. 

Summarizing, the procedure of calculating S{nB, np) consists of three steps. 
First, we pull out one bosonic (fermionic) brick from the left and right states. 
There is in general — 1 ( d{np) ) ways to do this. Second, we commute 
(anticommute) these two bricks, and replace the commutator (anticommutator) 
by a normally ordered operator. Third, we evaluate the matrix elements of this 
operator between states with a lower number of quanta (see section [3. 3p . 

However, before we move to the description of the evaluation of matrix 
elements of operators, we have to tackle the problem of orthonormalization of 
the basis states and of the evaluation of the commutators of composite bricks. 
We do this in the following two subsections. 

3.2.1 Linear independence and orthonormalization 

The recursive approach produces a set of Fock states in which some states may 
be contained in several copies. Therefore, one has to implement a mechanism 
to remove such redundancy. Because of the recursive structure, such mecha- 
nism will have to deal only with states with a given number of bosonic and 
fermionic quanta. Hence, its computational effort is small compared to what 
would it be should the whole Fock basis be considered. Moreover, the remain- 
ing states must be orthonormalized. We use a procedure which realizes these 
two tasks in one step. It is done by a numerical diagonalization of the matrix 
of scalar products S(nB,np). Subsequently, eigenvectors with corresponding 
nonzero eigenvalues are retained. Since the Fock states are not normalized, the 
normalized eigenvectors have to be multiplied by the inverse of the square root 
of their corresponding eigenvalues. We group such vectors in a matrix denoted 
by R{nB,np). Note that R{nB,np) is not to a square matrix in general. Then, 
we can write 

R^{nB,np)S{nB,np)R{nB,np) = l_D(„B,«f)xD(«B,nF), 

where lD{nB,nF)xD{nB .np) i^ the unity matrix which rank is equal to the size of 
the subspace of the Hilbert space with hb bosonic and np fermionic quanta. 

*Such an anticommutator is equal to a sum of operators involving only bosonic creation 
and annihilation operators. In general not all of these operators will be normally ordered 
operators. However, any such operators can be brought to a normally ordered form by appro- 
priately ordering the creation and annihilation operators, which they are composed of, using 
the commutation rules eqs.([T)l. Such ordering will produce additional operators which have 
to be taken into account and which can also be brought to a normally ordered. Thus, any 
anticommutator can always be written in a normally ordered form. 
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3.2.2 Automatic evaluation of commutators and anticommutators 

We can conclude from eqs. ([38]) and (p9|) that, in order to calculate a scalar 
product or a matrix element of some operator, the set of commutators and an- 
ticommutators of all composite bricks must be supplied. These commutators 
and anticommutators must be brought to a gauge invariant, normally ordered 
form and should be maximally reduced using the Cayley-Hamilton theorem. 
In general, there will appear some new single trace operators containing both 
creation and annihilation bosonic operators. Their commutators and anticom- 
mutators with all composite bricks should be evaluated as well and supplied to 
the algorithm. The number of such relations to be calculated grows rapidly, 
both, with increasing N, and increasing fermionic occupation number. Already, 
for the SU{3) gauge group one needs about a thousand of (anti)commutators. 
Therefore, a computer program was written to evaluate them. It uses standard 
relations H]), (U) and between the bosonic and fermionic creation and an- 
nihilation operators to move them among and within traces. As an example 
Eas. ([iD|) present one of the commutators needed for the calculations, 

[{a''a'iaa''a''aa), (a^a^)] = 1.5555{a^ a^) + 1. 3333 {a^ a'') {a^i a'' a) 
+ {a^a''a^){a''a''aa) + l.Ull{a'i a'' a''){a'i a) 
+ 0.l666(a1'at)(a'fa^at)(aa) -t- {a'' a^){a'' a'' aa^ a) 
-f 0.25(ata1')(atat)(a^aa) (40) 

We have found that it is more efficient to work with composite bricks than el- 
ementary ones. For example, if we have treated the operator {a' a' f^){a'^ p){a^ f'' f''), 
which is one of the composite bricks in the np = 4 sector for the SU(3) group, 
as a product of three operators, we would relate the desired matrix element 
with matrix elements from the np = 3 sector, which, in turn, would be related 
to some matrix elements in the np = 2 sector, and so on. However, one could 
treat it as a single operator, and jump to the bosonic sector in one step, consid- 
erably decreasing the number of commutation and anticommutation relations 
needed for such computations. On the other hand, the more elementary bricks 
will be contained in a composite brick, the more complicated the commutation 
and anticommutation relations will be. Hence, the usefulness of using more 
complex composite bricks is a question of balancing between the processor time 
consumed for evaluation of these relations, the memory needed to store them 
and the processor time gained by using more complex (anti)commutators. 

3.3 Matrix elements of gauge invariant operators 

In this section we describe the evaluation of matrix elements of operators, such 
as those appearing in the right hand side of eqs. ([55)1 and ([55]) . 

Let us denote a generic operator by 0(n^,n^). The arguments of O, n§ 
and Up have the following meaning 
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• n§ is the difference between the number of bosonic creation and annihi- 
lation operators contained in O, 

• n'p is the difference between the number of fermionic creation and anni- 
hilation operators in O. 

In general, n§ and can be any integers. Furthermore, we denote by 
0{ng,np)nB.nF the matrix element of O between basis states containing ub 
bosonic and n p fermionic quanta on the right hand side, and n'g = ng + rig 
bosonic and n'p — np + rip fermionic quanta on the left hand side, 

0(ng,ngU,„, = {{n's,n'p}\Oin'^,n'^)\{nB,nF}). (41) 

0{nB,np)nB,nF is a matrix of sizes D{nB,nF) x D{n'g,n'p), where the num- 
bers Dn^m denote the multiplicity of basis states with n bosonic and m fermionic 
quanta, as it was introduced in section [2.4.1l 

We first deal with some 'boundary' situations, and then we consider the 
generic case. 

3.3.1 Boundary cases 

Let #(a;)o denote the number of occurrences of the operator x in the operator 
O. The following observations can be exploited to simplify the computations: 

• if a matrix element of an operator for which #(/^)o > #(/)o is to be 
calculated then it is more convenient to evaluate its complex conjugate. 
Similarly, if we have to compute a matrix element of an operator for 
which = but #(a^)o > #(a)o we should rather compute 
its complex conjugate. 

• the matrix element of an operator O which has fermionic or bosonic an- 
nihilation operators acting on the Fock vacuum vanishes. 

• the matrix element of the bosonic elementary brick between states from 
the bosonic sector can be read off from the appropriate part of the matrix 
of scalar products. 

• the matrix element of an operator which is a product of two trace oper- 
ators can be calculated by inserting an identity operator between them, 
evaluating their matrix elements separately, and eventually multiplying 
and summing the partial results. 

3.3.2 Generic case 

In the generic case, we can assume that O is normally ordered (see footnote U] 
on pagell9p. If it is composed exclusively of creation operators, one can express 
it in terms of bosonic elementary bricks and use the appropriate boundary case, 
described above. 
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The strategy to evaluate a matrix element of O is to drag O over the op- 
erators constituting the right hand side state so that it annihilates the Fock 
vacuum. We start by pulling it through the composite fermionic brick, 

{0{n%,n$)ns,7i.p) ^ = {{n'B,n'p]\0{n%,n%)C\np,nF,p)\{nB,nF]) 

= {{n'g,n'p]\[0{n^,n^),C\np,nF,p)\\{nB -np,0}) 
+ {{ng,np]\C'^{np,nF,p)0{n%,n^)\{nB - ?^p,0}). 

(42) 

In order to move further, we substitute the relation for the (anti)commutator 
of 0(n§,n^) and C^(p, For each operator appearing in this relation 

we evaluate its matrix element, first checking whether conditions for any of the 
special cases are met. This task should be easier, since these matrix elements 
must be evaluated between states with smaller number of bosonic and fermionic 
quanta. The second term in eq. (1421) . can be calculated by inserting an identity 
operator between the operators C'^ [up^np^p) and 0{ng,np). Again, this com- 
putation should be easier. On one hand, the matrix element of the elementary 
brick C'' {np,np,p) should be known from the matrix of scalar products. On 
the other hand, the matrix element of 0{n'^,np) involves states with smaller 
number of bosonic and fermionic quanta. 

The purely bosonic case can be treated analogously. We have 

(0(ng, 0)„,,o) = {{n'B,0} I [0(ng, 0), 0)] \{nB - p, O}) 

+ ({n^,0}|C(p,0)O(ng,0)|{nB -p,0}), (43) 

In principle, we can proceed with those relations until O hits the Fock vac- 
uum. Since O is a normally ordered operator, such matrix element vanishes by 
definition. Therefore, collecting all the intermediate results, we should be able 
to evaluate the desired matrix element. Nevertheless, one has to remember that 
states in Ij/is, O}) are in general not orthonormal, so one has to implement into 
relations (|42|l and (|43|) the orthonormalization procedure. 

3.3.3 Recurrence relations 

Since the recursively constructed set of basis states can contain degenerate 
states, the matrices of operators, calculated as described above, will also con- 
tain such redundant matrix elements. In order to get rid of them one has to use 
the R{nB,np) matrix. Incorporating this matrix in the relations for the ma- 
trix element of any operator, enables us to formulate the complete and correct 
recurrence relations. We have two recurrence relations: 

• the expression of the matrix element in the fermionic sectors in terms of 
matrix elements in the bosonic sector, 
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{{Ti'B,n'p}\0{n'^,n'^)\{nB,TiF}) = 

(({"'b, [0(ng, nO), nF,p)]\{nB - n^, O}) 

+ {{nB,np}\C^np,nF,p)0{n%,np)\{nB -~ np,0})j • R{nB,nF), 

(44) 

and 

• the expression of the matrix element in the bosonic sector with bigger 
number of bosonic quanta in terms of matrix elements with smaller number 
of bosonic quanta, 

({n's, 0}|O(ng, 0)|{nB, 0}> - (({n^, 0} | [0(ng, 0), Ct(p, 0)] [{ns - p, O}) 

+ {{n'B,0]\C\p,0)O[n%M{nB-p,0]))-R{nB,Q). (45) 

Thus, we have expressed the desired matrix element in terms of matrix elements 
of operators between states will lower number of quanta, which should have been 
already evaluated during some previous calculations. 

With these recurrence relations the presentation of the whole algorithm is 
complete. One can use them to evaluate the matrix of scalar products in the 
sectors which have not been considered so far in the calculations. This done, 
the matrix elements of operators, needed for the calculations in sectors with yet 
bigger number of quanta, can be computed. In this way one can proceed until 
the cut-off Ncut is reached. 

Note that the above algorithm is very universal. In principle it can be used to 
systems defined in space of any dimensionality. Particularly, if the Hamiltonian 
is invariant under a SO{d) symmetry our method can be generalized in order 
to calculate the spectra in the channels with given angular momentum [12j . 
Moreover, it is applicable to systems with bosonic and fermionic polynomial 
interactions as well as to systems with discrete or continuum spectrum. 

In the next section we describe the results for supersymmetric model with 
discrete spectrum obtained with our approach. 

4 Applications 

As an application of the above algorithm we present results for a supersymmetric 
system given by the Hamiltonian constructed in 12.21 

H^i,p'+g^{trx^-^{tTx^f)^2gir{x[f\f]), (46) 

where the gauge symmetry group is chosen to be the SU{'i) group. 

The main motivation for studying this Hamiltonian is that, being a system 
with a discrete spectrum, it is an good test-ground for our numerical method. 
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Moreover, the bosonic part of its potential is similar to the potential of SYMQM 
with the tensor fate replaced by date [12 • Therefore, one hopes that some 
analytic approaches based on the numerical results from such simple model can 
be tested for a future application to the more dimensional SYMQM systems. 
On the other hand, eq. (B5|) is interesting by itself. Being an supersymmetric 
anharmonic oscillator it contains a nontrivial fermionic interaction. 

In the following we start by analyzing the convergence of the eigenvalues, 
then present the spectra calculated numerically. Subsequently, we briefly de- 
scribe their symmetries, such as supersymmetry and scaling symmetry. Even- 
tually, we calculate numerically the Witten index for this model. 

4.1 Numerical spectra 

One of the advantages of the cut Fock space approach as a numerical method is 
that it enables one to judge on the reliability of the results and estimate their 
errors. To this end, the convergence of eigenenergies or eigentates with increas- 
ing cut-off must be investigated. Table [7] contains the energies of the lowest 
eigenstate in fermionic sector with np = 0, ... ,4 obtained for different N^ut- 
The results in sectors with higher fermionic occupation number are related to 



Neut 


np = 


Up = 1 


np = 2 


1 


2.833333333333333 


3.750000000000124 


3.75000000000008 


5 


2.805137759654418 


2.817654396966426 


2.41010649311797 


10 


2.804878933491876 


2.804943385906189 


2.38393952020263 


15 


2.804877899477374 


2.804878578502977 


2.38379874405844 


20 


2.804877857980324 


2.804877869314702 


2.38379576457689 


25 


2.804877857812559 


2.804877857890121 


2.38379573799721 


30 


2.804877857802534 


2.804877857804384 


2.38379573773261 


35 


2.804877857802529 


2.804877857803605 


2.38379573772474 


40 


2.804877857802507 


2.804877857803596 


2.38379573772458 


Ncut 


np ~ 3 


np = 4 




1 


1.009109012532963 


3.750000000000082 




5 


0.017808308382480 


2.001903558629864 




10 


0.000102896003680 


1.978068350121234 




15 


0.000002580129746 


1.977963562445323 




20 


0.000000018188518 


1.977960963230507 




25 


0.000000000380532 


1.977960939859982 




30 


0.000000000013940 


1.977960939644200 




35 


0.000000000004579 


1.977960939638051 




40 


0.000000000004250 


1.977960939637698 





Table 7: The dependence on the cut-off N^ut of the lowest eigenenergy in several 
fermionic sectors. Note the exact degeneracy of the levels in the np = and 
np — I sectors due to supersymmetry and the appearance of the supersymmetric 
vacuum in the np = 3 sector. 
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the ones presented in the table through the particle-hole symmetry. A conver- 
gence to more than 10 digits is achieved with N^ut — 40, which corresponds to 
a Fock basis of about 150 states in the bosonic sector and about 1300 states in 
the = 4 sector. Therefore, one can safely use the cut-off Ncut = 40 in order 
to evaluate the lowest eigenenergies. The uncertainties of those eigenenergies, 
defined as the difference of the outcomes for consecutive cut-offs, are negligible. 
Figure [T] and figures [5] show the dependence of 15 lowest eigenenergies on the 
cut-off in different fermionic sectors. All the figures were made for g = 1.0. 




6 - z 

4 - 

2 I ^ ^ ^ , ^ , ^ 

5 10 15 20 25 30 35 40 

Ncut 

Figure 1: Convergence of the 15 lowest eigenenergies in the bosonic sector with 
increasing cut-off. An exponential-like convergence is seen, which is a charac- 
teristic feature of bound states. 

We can clearly see a fast, exponential-like, convergence, which was shown to 
be a characteristic feature of bound states [2] [15]. Hence, one can conclude 
that indeed all the spectra are discrete. This should be contrasted with the 
results for the SYMQM systems. The latter have potentials with flat directions 
which induce continuum spectra in some sectors [H][I2][I1]. Table [S] contains the 
eigenenergies of lowest states in all 9 fermionic sectors calculated with g = 1.0 
and a cut-off Ncut = 40. We comment these results in more details in the 
following subsection. 
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4.2 Symmetries 
4.2.1 Supersymmetry 

Supersymmetry can be seen in table |8] as an exact degeneracy of the converged 
eigenenergies in neighboring fermionic sectors. One notices a non-degenerate 
vacuum state in sector with np = 3 and its image through the particle-hole 
symmetry in the sector with np = 5 signaling an unbroken supersymmetry. 
The number of supersymmetric vacua was discussed in [27J in the case of free 
SU{N) model. Accordingly, there are four zero-energy states in sectors with 
np — 0,3,5,8. Hence, two of these states disappear when the interaction is 
turned on. The particle-hole symmetry is responsible for a double degeneracy 
of eigenenergies in the sector with np — A. This sector must contain as many 
states which are parts of supermultiplets formed with states from the sector 
with nF = 3 as there are states which form supermultiplets with states from 
the sector with np = 5. Moreover, both sets of supermultiplets must have the 
same spectra. This is a specific feature of models with SU(N) gauge group with 
N odd. 

Table [7] reveals also signatures of supersymmetry. One one hand, the lowest 
states in sectors with np = and np = 1 are exactly degenerate. On the 
other hand, the presence of a non-degenerate supersymmetric vacuum state in 
the sector with np = 3 and in the sector with np — 5 is a proof that the 
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np = 2 


np = 3 


np = A 


np = 5 
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np = 7 


np = 8 
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1.97796^ 
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4.98444 
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5.75469 
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5.75469 




5.75469 


5.75469 








5.75469 


5.75469 
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5.75469^ 


5.75469 
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5.85572^ 
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6.21387 


6.21387^ 


6.21387 












6.48704 


6.48704 




6.48704 


6.48704 












6.82769 


6.82769^ 


6.82769 














6.88508 


6.88508^ 


6.88508 












6.89869 


6.89869 




6.89869 


6.89869 






7.49223 


7.49223 












7.49223 


7.49223 






7.66421 


7.66421 




7.66421 


7.66421 












7.75149 


7.75149^ 


7.75149 














7.81189 


7.81189^ 


7.81189 












7.82955 


7.82955 




7.82955 


7.82955 










8.11747 


8.11747 




8.11747 


8.11747 






8.39012 


8.39012 












8.39012 


8.39012 




8.39012 


8.39012 








8.39012 


8.39012 










8.39012 


8.39012^ 


8.39012 









Table 8: Eigenenergies of few lowest eigenstates with g 1.0. The states in 
the np — A sector are double degenerate which follows from the particle- hole 
symmetry. Hence the ^ notation. An exact, supersymmetric pairing can be 
observed among states in adjacent fermionic sectors. 



supersymmetry is unbroken in this model. 
4.2.2 Scaling symmetry 

The scaling property of the quantum anharmonic oscillator was first noted by 
Symanzik, and elaborated by Simon |28) . If we consider the transformations 

x^Xx, P^\p, / -> /, P -> /^ (47) 
then the Hamiltonian is rescaled as 

H^H=^Up'+ g^X^ (tr x' - l(tr x^f) + 2g X' tr(a;[/t, /]) ] . (48) 
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(49) 



Since, the transformations ea. (l47p can be unitarily implemented, both Hamiho- 
nians in ea. (|15|) have identical eigenvalues. Therefore, it is sufhcient to calculate 
the spectrum at g = 1.0. 

Figure |3] shows the dependence of six lowest eigenenergies from the bosonic 
sector on the coupling constant. The numerical results are compared with the 
prediction of eg. (115)) . For large g the agreement of both should be noted, whereas 
the discrepancies for small g are due to finite cut-off effects as discussed below. 

4.3 Critical slowing down 

A critical slowing down can be observed on figure [3] in the vicinity oi g — 0. For 
a small value of the coupling constant a much higher cut-off is needed in order 
to obtain converged results, ft is because aX g — the spectrum is free and the 
eigenenergies corresponding to nonlocalized states calculated by our algorithm 
do not converge. Rather, they fall off to zero with increasing cut-off in a power 
like manner, a behavior resulting from approximating a plane wave by a finite set 
of localized harmonic oscillator eigenstates. Hence, the eigenenergies calculated 
for small coupling constant at small cutoff cannot follow the curve eq. (j49|) of 
exactly converged energies. Nevertheless, with increasing coupling constant the 
eigenenergies of bound states converge more and more rapidly. Those energies 
that have already converged with the cut-off agree with the analytic prediction 
of eq. (|1^ . The complete discussion of such critical behavior is out of scope of 
the present work and will be carried out elsewhere. 

4.4 Perturbative expansion 

In a cut Fock space the states corresponding to the free states in the continuum 
are normalizable. Moreover, their exact analytic form can be obtained for a 
finite cut-off |30| |31j. Therefore, at fixed finite cut-off N^^t one can perform a 
perturbative expansion in the coupling constant g. Using well-known formulae 
for the perturbative corrections to the energy [32] one can obtain the approx- 
imate dependence of the eigenenergies on the coupling constant. We present 
here the results for the lowest eigenstate from the bosonic sector calculated for 
the cut-off Ncut = 20, 



E = Eo + g^Vs, 



'0 




+ 



+ 9' 




^Eo,E' 



) 



{Eo - E'Y 



(50) 
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0.1 0.2 0.3 0.4 0.5 0.6 

9 

Figure 3: Dependence of the first six eigenenergies on the couphng constant g in 
the bosonic sector. Crosses denote the numerical results obtained with cut-off 
Ncut — 40, whereas the dashed lines correspond to the prediction of eq. (|49l) 
with H(\) evaluated numerically. The finite cut-off effects can be seen for small 
values of the coupling constant. 



where Ve^e' is the matrix element of the potential between states with energies E 
and E' respectively. The comparison of this formula with the numerical results 
and the prediction of ea. (|l^ is presented on figure U] It is remarkable that the 
perturbative expansion attains the values of the coupling constant where the 
numerical results have converged for cut-off Ncut and agree with eq. (|49|) . This 
gives us hope that with an improved perturbative expansion the whole spectrum 
may be obtained analytically. This issue is being investigated in more details. 

4.5 Witten's index 

The Witten index is defined as 

/^(r)= E (51) 

6Gbosoiiic states /Gfermionic states 

It is a commonly used quantity to study supersymmetry in quantum mechanics. 
In our set up, the sums over the bosonic and fermionic states are finite due to 
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Figure 4: Comparison of the dependence of the energy of the lowest bosonic 
eigenstate on the couphng constant g. Numerical data (crosses) are plot to- 
gether with the prediction of eq.(|49p (solid line), and the perturbative expansion 
eq. ([50|) fdashed lines) 



the cut-off. For a given Ncut one can just plug the eigenenergies obtained from 
each fermionic sector into eq. ()51|) . Our numerical results are shown on figure 
[SI where the dependence of Iw{T) on the euclidean time T is presented. One 
notices a rapid convergence of Iw{T) to the value —2. This confirms the fact 
that the model has two supersymmetric vacua. 

The collapse of the Witten index at T = to zero can be explained in the 
following way. At T = the Iw{T) is just the difference in the numbers of 
bosonic and fermionic states. However, the cut Fock basis has an equal number 
of bosonic and fermionic states (see subsection 12. 5p . Thus, the value of the 
Witten index at T = is zero for any cut-off. This is no longer true at T 7^ 0. 
Especially, for T — >■ cxd, contributions from the states with nonzero energies 
cancel, and Iw{T) counts the numbers of supersymmetric vacua, equal —2 in 
the present case. 
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Figure 5: Witten index. Different curves correspond to increasing cut-offs: 5, 
10, 15, ... , 40. 



5 Conclusions 

In this paper we have described in a very detailed manner a recursive algorithm 
for evaluation of matrix elements of any gauge invariant operator in a Fock basis 
of Hilbert space. It can be applied to systems with any SU{N) gauge group 
and allows the evaluation of spectra in all fermionic sectors. We demonstrated 
the calculations on an example of anharmonic oscillator with supersymmetric 
interactions for the SU{3) group. 

We started by describing the idea of the numerical approach to quantum me- 
chanics in the Hamiltonian formulation using the cut Fock space method. Then, 
we presented the construction of gauge-invariant Fock basis, and particularly we 
introduced the concepts of elementary bosonic bricks and composite fermionic 
bricks. Such approach provided us a systematic, recursive description of the 
Fock states with increasing number of quanta. We discussed the properties and 
symmetries of such basis. Next, we concentrated on the numerical algorithm. 
We described the calculation of the matrix of scalar products, emphasizing the 
main ideas of the recursion relations. Then, other parts of the algorithm were 
presented: the procedure which removes redundant basis vectors and orthonor- 
malizes the remaining ones as well as the program which automatically calcu- 
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lates the commutators and anticommutators of given operators. Finally, the 
expressions for the evaluation of matrix elements of any operator were outlined, 
and eventually, the full recursion relations were presented. In the third part of 
this article we applied our algorithm to a supersymmetric system with SU{3) 
gauge group and a discrete spectrum. We used it as a particularly well suited 
test-ground for our approach. We calculated the eigenenergies in all 9 fermionic 
sectors and discussed their symmetries. Eventually, we were also able to obtain 
the Witten index for this system. 

The main advantage of this algorithm is that it treats bosons and fermions 
on an equal footing, and thus, enables calculations in any fermionic sector of the 
Hilbert space. This should be contrasted with the sign problems encountered in 
lattice field theories. As a result an exact supersymmetric degeneracy can be ob- 
tained even for finite cut-off. Moreover, the approach can be applied to systems 
with discrete and continuous spectra as well as possessing any kind of gauge 
symmetry. Particularly, the supersymmetric anharmonic oscillator presented in 
this article with gauge groups with > 3 is currently investigated. Similarly, 
the Z? = 2, supersymmetric Yang-Mills quantum mechanics with several gauge 
groups, such as SU{3), SU{A) and SU{5) are studied analytically and numeri- 
cally [30] [31] . The flexibility of the algorithm enables also an generalization to 
higher dimensions, with the ultimate D = 10, SYMQM case in mind. Results 
of the work is this direction are promising. 

Acknowledgments 

The Author would like to thank for discussions with prof. J. Wosiek on the 
subject of this paper and his careful reading of the manuscript. 

References 

[1] T. Banks, W. Fischler, S. Shenker, L. Susskind, 'M-theory as a matrix 
model: a conjecture', Phys. Rev. D 55 (1997) 6189 

[2] J. Hoppe, 'Quantum theory of a massless relativistic surface and a two 
dimensional bound state problem', PhD thesis MIT, 1982, unpublished 
(scanned version avaible at |http://www.aei-potsdam.mpg.de/| hoppe) 

[3] M. Liischer, 'Some analytic results concerning the mass spectrum of Yang- 
Mills gauge thoeries on a torus', Nucl. Phys. B 219 (1983) 233-261 

[4] M. Liischer, G. Miinster, 'Weak-coupling expansion of the low-lying energy 
values in the SU{2) gauge theory on a torus', Nucl. Phys. B 232 (1984) 
445-472 

[5] P. van Baal, 'Gauge theory in a finite volume'. Acta Phys. Pol. B 20 (1989) 
295 



32 



[6] K. N. Anagnostopoulos, M. Hanada, J. Nishimura, S. Takcuchi, 'Monte 
Carlo Studies of Supersymmetric Matrix Quantum Mechanics with Sixteen 
Supercharges at Finite Temperature', Phys. Rev. Lett. 100 (2008) 021601 

[7] M. Hanada, A. Miwa, J. Nishimura, S. Takeuchi, 'Schwarzschild radius 
from Monte Carlo calculation of the Wilson loop in supersymmetric matrix 
quantum mechanics', Phys. Rev. Lett. 102 (2009) 181602 

[8] J. Wosick, 'Spectra of supersymmetric Yang-Mills quantum mechanics', 
Nucl. Phys. B 644 (2002) 85-112 

[9] J. Wosiek, 'Supersymmetric Yang-Mills quantum mechanics in various di- 
mensions', PoS LAT 2005 (2006) 273 

[10] J. Wosiek, 'Spectra of supersymmetric Yang-Mills quantum mechanics', 
Nucl. Phys. B 644 (2002) 85-112 

[11] M. Campostrini, J. Wosick, 'Exact Witten index in D=2 supersymmetric 
Yang-Mills quantum mechanics', Phys. Lett. B 550 (2002) 121-127 

[12] M. Campostrini, J. Wosiek, 'High precision study of the structure of D=4 
supersymmetric Yang-Mills quantum mechanics', Nucl. Phys. B 703 (2004) 
454-498 

[13] P. Korcyl, 'Classical trajectories and quantum supersymmetry', Phys. Rev. 
D 74 (2006) 115012 

[14] M. Trzetrzelewski, J. Wosick, 'Quantum systems in a cut Fock space'. Acta 
Phys. Polon. B 35 (2004) 1615-1624 

[15] M. Trzetrzelewski, 'Quantum mechanics in a cut Fock space', Acta Phys. 
Polon. B 35 (2004) 2393-2416 

[16] P. Korcyl, 'Eigenvalues and eigenvectors of the d dimensional Laplace- 
Beltrami operator in a cut Fock basis', in preparation 

[17] F. Cooper, A. Khare, U. Sukhatme, 'Supersymmetry and Quantum Me- 
chanics', Phys.Rept. 251 (1995) 267-385 

[18] M. Trzetrzelewski, 'Reduction of SU{N) loop tensors to trees', J. Math. 
Phys. 46 (2005) 103512 

[19] M. Trzetrzelewski, 'Supersymmetric Yang-Mills quantum mechanics with 
arbitrary number of colors', Ph.D. thesis, Jagiellonian University 

[20] M. Trzetrzelewski, 'Number of gauge singlets in supersymmetric Yang-Mills 
quantum mechanics', Phys. Rev. D 76 (2007) 085012 

[21] Hamermesh, 'Group theory and its application to physical problems', 
Addison- Wesely, Reading Mass, 1962 



33 



[22] M. Byrd, 'The geometry of SU{3)\ LANL ePrint pysics/9708015| (1997) 



[23] T. Tilma, E.C.G. Sudarshan, 'Generalized Euler Angle Parametrization for 
SU{N)\ J. Phys. A: Math. Gen. 35 (2002) 10467-10501 

[24] H. Weyl, 'The Theory of Groups and Quantum Mechanics', Methuen & 
Co. Ltd., 1931 

[25] S. Samuel, 'Solutions of extended supersymmetric matrix models for arbi- 
trary gauge groups', Phys. Lett B 411 (1997) 268-273 

[26] M. Claudson, M.B. Halpern, 'Ground state wave functions', Nucl. Phys. B 
250 (1985) 689-715 

[27] M. Trzetrzelewski, 'Supersymmetry and Lie groups', J. Math. Phys. 48 
(2007) 083508 

[28] B. Simon, 'Coupling constant analycity for the anharmonic oscillator', Ann. 
of Phys. 58 (1970) 76-136 

[29] P. Korcyl, 'Study of a supersymmetric quantum system with discrete spec- 
trum for different gauge groups', in preparation 

[30] P. Korcyl, 'Exact solutions to D — 2, Supersymmetric Yang-Mills Quantum 
Mechanics with SU{3) gauge group'. Acta Phys. Pol. B 2 (2009) 623, arXiv: 
0911.2152 

[31] P. Korcyl, 'Exact solutions to D — 2 Supersymmetric Yang-Mills Quantum 
Mechanics using a cut Fock space method', in preparation 

[32] L. D. Landau, E. M. Lifschitz, 'Quantum mechanics: non-relativistic the- 
ory', 2cond edition 

[33] E. Witten, 'Dynamical breaking of supersymmetry', Nucl. Phys. B 185 

(1981) 513-554 

[34] E. Witten, 'Constraints on supersymmetry breaking', Nucl. Phys. B 202 

(1982) 253-316 



34 



